General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



OST t NX 



RN-78 

ANNUAL PROGRESS REPORT 


NASA GRANT NGL 22-009~b48 
"Flight Control Systems Research" 

1 September 1974 to 31 August 1975 


Principal Investig^'tor : 


Professor H. Philip Whitaker 


/ 


I 


MASSACHUSETTS INSTITUTE OF TECHNOLOGY 
Department of Aeronautics and Astronautics 
Measurement Systems Laboratory 
Cambridge, Massachusetts 02139 


T.ch.) ,] 

• ^ CSCT T 1 K 


ci/n 


V7«5-3111» 


”ncl as 
^5325 


■c. 


AC.’<NOWLEDGEMENT 


This report was prepared under DSR Project 73598 sponsored 
by the NASA-Edwards Flight Research Center, through National 
Aeronautics and Space Administration grant number NGL-22-009-548. 

The publication of this report does not constitute approval 
by the NASA of the findings or conclusions contained herein. It 
is published only for the exchange and stimulation of ideas. 


Annual Progress Report 
NASA Grant NGL 22-009-548 

Prepared by 
H. P. Whitaker 

During this reporting period under the step funding provision 
of the Grant, manpower was reduced and reached a level of 25% of 
the principal investigator in February 1975. There has been a 
corresponding reduction in output. The main emphasis of the 
activity has been placed upon evaluating and testing the most re- 
cent method of calculation of the performance index, application 
to various example designs, and progreunming for the discrete data 
control system case. The specific items are summarized in the 
paragraphs that follow. The detail descriptions are to be con- 
tained in a Contractor Report, the writing of which will form a 
major part of the activity under the Grant in the coming year. 

1. Computation of the performance index* 

Due to the inaccuracies that were introduced by digital word 
length round-off and by the matrix inversion routines, an alter- 
native method of computing the performance index was developed. 
This involves the direct integration of the performance index 
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and makes use of the partial fraction expansion of the Laplace 
transform of the excitation input to the error state. This 
technique requires accurate knowledge of the poles and zeroes 
of that transform so that accurate values oJ residues can be 
obtained. . This in turn made it necessary to spend effort modi- 
fying that portion of the program, and finally a program package 
for finding the eigenvalues of matrices prepared at the Argonno 
Laboratory was adopted. In tho examples investigated to this 
date, the performance index calculation has exhibited acceptable 
accuracy. A brief development of the new computational techni- 
que is presented in Appendix A. 

2. Corresponding changes were made to the computer program for 
the discrete design case. That progrsun has not been tried out 
to the extent that the program for the continous systems has. 

There remain several areas that need to be modified. It is not 
anticipated that these involve major programming efforts, but 
rather involve changes that are needed to expedite the manner in 
which data is presented in the input description. Since both con- 
tinuous and discrete sub-systems need to receive input data 
description, the input requirements should be simple and easy to 
implement to avoid input data errors. 

The discrete program has been checked out using the same 
configuration that was evolved for an analog C* system. Figure 1 
presents a mathematical block diagram for the system. The air- 
plane and elevator servo are analog, and the rest of the control 
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Kathen«tical Block Diagram of an Aircraft C* Control Syst 
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2. Step C* Comnand Response 


system configuration is provided by the digital section. The 
C* error signal is formed and fed to the elevator. Inner loop 
feedbacks of pitch rate, h, and normal acceleration, a^, are 
provided to give a closed loop static sensitivity of unity. 

Figure 2 presents the step function C* response. The optimiza- 
tion of the continuous system indicated that the pitch rate feed- 
back sensitivity ratio was insensitive, and for this preliminary 
check of the digital program that gain was kept constant. The 
variable parameters were chosen to be the forward loop static 

sensitivity, the acceleration feedback, and 

c IUJ6 j c ia2 ,0 j 

the integrator sensitivity, S^. The model -esponse shown in 
Figure 2 is that of the discrete model which differed from the 
continuous model by less than 5% due to the difference between 
an analog and digital representation. The model corresponds to 
a second order lag with - 9.5 rad/sec and X • 0.5. Table 1 
summarizes the parameter values that result- ..d. The optimization 
was terminated 



Tab le 1 


Parameter 

Analog System 

Discrete System 


2.02 

1.84 


o^ 
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1 

-1.13 

sec”^ 

0.504 

0.504 


when the parameters become insensitive. It may be that even 
closer agreement would be found at the mathematical optimum. 
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Figure 2 indicates that such variation is insignificant however 
in a practical design. 

3. Search algorithm: The algorithm for finding the minimum of 

the performance index is a straightforward gradient search. Some 
control of step size is incorporated, and effort was expended in 

f 

improving the convergence characteristics as more design examples 
have been investigated. No algorithm seems to be optimum for 
all cases, but the present version appears to be a practical com- 
promise. Appendix B presents a summary of the algorithm. 

4. Lateral control system studies: As an example of a more com- 

plicated control system configuration, an investigation of a 
rudder coordination system for the F-8 was undertaken. One can 
select several different indications of desired coordination: 

(a) the minimization of the excitation of the lateral oscillation 
for aileron Inputs, (b) minimizing sideslip by causing the yaw 
stability axis component of angular velocity to be proportional 
to roll angle during the entry into a constant bank angle turn, 
or (c) minimizing the average squared sideslip during turn entries. 
Cases (a) and (b) involve different model choices, and case (c) 
involves a slightly different performance index. Since the pre- 
sently used version of the performance index computation evaluates 
an expression of the form 

(PI) - /[i(t)]*dt 
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where the Laplace tranaform of l(t) is known in pole-zero form, 
the use of 

(PI)g - /lB(t)]»dt 

is a minor modification requiring only that B rather than i be 
taken as the output quantity. When case (b) was examined, the 
closed-loop zeroes of the yaw rate transfer function apparently 
cause a large initial transient effect that leads to numerically 
high values of the minimum performance index. This seems to be 
unsatisfactory and needs to be investigated more thoroughly. 
Otherwise all three examples gave satisfactory results. 

Planning 

The next year of this Grant will be devoted to writing the 
final report and to more thorough examination of the discrete 
design case. During the latter, emphasis will be placed upon the 
effect of sampling frequency, and preliminary work on this has 
already begun. 


5 


Appendix A. Computation of the Performance Index 
Reference: 

Palsson, T. , Parameter Uncertainties in Control System 
Design , Measurement Systems Laboratory Report TE-46, 
May 1971. 


Section 3.5 of the referenced report shows that the Model 
Performance Index is directly equivalent to 

(PI) - /Ii(t)]*dt 


where i(t) is the excitation input to the dynamic model of the 
error states. Referring to figure A-1, the error state can be 
obtained in two ways, giving the two expressions 
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Solving for i 




n+k 

i (p) « E R./(p-p. ) 

i-1 ^ ^ 


where p^ are the (n + k) poles obtained from the system poles 
and the model zeroes; and R^ are the corresponding residues. 
Then 

n+k ^ 

i(t) = Z R.ePi' 
i»l ^ 
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Real (p^) < 0 


Thus by obtaining the partial fraction expansion of i (p) 
one can evaluate the perfonnance index explicitly. This also 
requires having the eigenvalues of the augmented system. 
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Appendix B. Currently Used Search Algorithm for the Parameter 
Optimization Program 7 August 1975 

This memo summarizes the algorithm used to search for the 
minimum of the performance index as of the modification made up 
to 22 July 1975. 

1. The changes made to the parameters at the end of an 
iteration are proportional to the gradient of the performance 
index. The proportionality factor, SGRD, is initially taken to 
be 10.0. The predicted change to the performance index is cal- 
culated using the parameter increments, and a correction factor, 
CORK, is calculated which is the ratio of the desired fractional 
change in the performance index to the predicted fractional 
change. In effect the proportionality factor is then corrected 
by multiplying by CORR. If the first order p*‘ediction were 
accurate, the actual fractional change in the performance index 
on the next iteration would be equal to the desired change. 

2. If a parameter has reached a point at which its gradient 
changes sign, its contribution to the reduction of the perfor- 
mance index is reduced by the step factor, STF, wn^ ;h is generally 
equal to or less than one. The gradient step change is multiplied 
by STF. When a parameter's gradient changes sign, its step fac- 
tor is reduced to 0.6 of its previous value. Thus if a parameter 
is alternately stepping across its optimum value, its permif-ed 
change decreases rapidly. To allow for the case for which the 
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gradient changes sign as a result of other parameters reaching a 
more favorable point In parameter space, the step factors are 
allowed to recover by Increasing them by 15% when success jvi.> 
Iterations result In the same sign of the gradient. The Increase 
is not applied If the step factor Is 1.0 or greater. 

3. When the gradients are small, the predicted change In 
parameters needed to obtain a desired change In the performance 
index may be large. The correction factor Is further modified 
so that no par 2 uneter will change by more than 200% In any given 
Iteration. 

4. After the resulting parameter changes have been cal- 
culi' id. the parameters are updated, and any parameter con- 
straints applied. 

5. The search algorithm also provides for automatic adjust- 
ment of the called for change in the performance Index on any 
Iteration. This takes place by modifying the parameter DELTA 
whose inltidl value Is read as part of the input data. DELTA is 
the desired fractional change In (PI) in a given iteration. Since 
(PI) seldom varies linearly with a parameter value, the linear 
prediction is good only far from the optimum point where the slope 
changes slowly. By modifying DELTA one can slow the parameter 
change as the minimum is approached and cut down the tendency to 
overshoot the minimum. The success ratio, SR, is defined as the 
ratio of the actual change in (PI) to the predicted change. 
Typically a success ratio of 0.5 is good. The program increments 
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DELTA by 120% whe.iever SR is outside the range 0.4 to 0.8. This 
is done in the OPTl subroutine. 

6. When the mlniinuin is approached, large changes in (PI) 

are not possible. If one has specified a value other than zero for 
the reference (PI) value, DELTA is further changed so that the 
called for change in (PI) would not be greater than one that 
changes (PI) co the reference value. This again prevents excessive 
parameter changes near the optimum point. It may not be necessary 
however. 

7. Further logic affecting the search is provided in the OPTl 
cubroutine. If the parameter change causes the system to become 
unstable, all parameter changes are decreased to 40% and the itera- 
tion re-started. If thr parameter changes cause the performance 
index to increase over the value on the previous iteration, the 
changes are reduced by 70%. If this still results in Increased 
(PI), the maximum parameter change is reduced to 5%. If this still 
results in passing the minimum, the search terminates since a 5% 
change in the parameters should be within practical tolerance 
bounds on the design parameters. One is cautioned however that if 
one parameter (or more perhaps) is very insensitive in this region 
of parameter space at which the search terminates, a large change 
in the insensi tive parameter may ta)te the process to a. region at 
which further optimization of the other parameter can result in a 
significantly lower value of the (PI) . The current progreim leaves 
it to the engineer to use his judgement to restart the parameter 
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if he feels further optimization is necessary rather than letting 
the gradient search proceed in very small increments using ex- 
cessive computation time. 

8. Various stopping conditions are used in the OPTl sub- 
routine. One can specify (PIREF) in the input data as a value of 
(PI) that is acceptadsly small and stop the search wht ' it is 
reached. If the change in (PI) is les£ than 5% on three neccessive 
iterations I the search terminates on the assumption that the DELTA 
modifications are indicating that the minimum point has been 
reached. It also terminates on reaching the specified maximum 
number of iterations. This can be used to calculate (PI) for 
given parameters values by inputting ITMAX * 0. If the normalized 
slope is less than 10% of the specified reference (PI) , it would 
ta)ce a ten- fold parameter change to make a change in (PI) equal to 
the reference (PI) . That is an indication of acceptably small 
gradients and hence that the minimum has been found, so the search 
is terminated. 
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